Transport properties of a multichannel Kondo dot in a magnetic field 
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We study the nonequilibrium transport through a multichannel Kondo quantum dot in the pres- 
ence of a magnetic field. We use the exact solution of the two-loop renormalization group equation 
to derive analytical results for the g factor, the spin relaxation rates, the magnetization, and the 
differential conductance. We show that the finite magnetization leads to a coupling between the 
conduction channels which manifests itself in additional features in the differential conductance. 
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Introduction. — The study of a localized spin coupled 
via an antiferromagnetic exchange interaction J to K in- 
dependent electronic reservoirs has a long history in con- 
densed matter physics. In the simplest case of if = 1 
the electron spins completely screen the local spin at low 
energies and thus lead to a Fermi liquid. In a renormal- 
ization group (RG) analysis this situation is characterized 
by the divergence of the renormalized exchange coupling 
J (A) at the Kondo scale A = Tk. The situation is com- 
pletely changed^ if the spin is coupled to more than one 
screening channel {K > 1). Then the renormalized ex- 
change co upli ng stays finite and flows to a non-trivial 
fixed point^ J* ^ 1/iC at low energies, which mani- 
fests itself in unusual non-Fermi liquid behavior like a 
non-integer "ground-state degene racy " or characteristic 
power laws in various observables.'^ 

The recent developments in the ability to engineer 
devices on the nanoscale lead to the experimental 
realizatiorP of two-channel Kondo physics in a quantum 
dot set-up.'^ In this set-up it was possible to measure 
the differential conductance and observe universal scal- 
ing and square-root behavior which are characteristic for 
the two-channel Kondo effect. This triggered theoretical 
studietP^fi^ of the transport properties of multichannel 
systems using conformal field theory as well as numerical 
and perturbative RG methods. The latter uses a pertur- 
bative expansion in the renormalized exchange coupling 
which is well-controlled provided K ^ \. Specifically, 
Mitra and Rosch- calculated the differential conductance, 
the splitting of the Kondo resonance in the T matrix, and 
the current-induced decoherence in the absence of a mag- 
netic field. Recently the spin dynamics was studiecP^ in 
the absence of a bias voltage and shown to possess pure 
power-law decay with an exponent g = A/ K. 

In this Brief Report we extend the analysis of Mitra 
and RosclP to include an external magnetic field. We 
perform a real-time RG (RTRG) analysis^-'^ to derive the 
renormalized magnetic field, the spin relaxation rates, 
the magnetization of the quantum dot, and the current. 
In particular, we focus on inelastic cotunneling processes 
which lead to characteristic features in the differential 
conductance whenever one of the applied bias voltages 
ViS equals the value of the renormalized magnetic field. 
We further show that the finite magnetization on the 



dot leads to a coupling of the conduction channels which 
results in additional features in the differential conduc- 
tance. 

Model. — We consider a quantum dot possessing a spin- 
1 /2 degree of freedom, which is coupled via an exchange 
interaction J to K independent electronic reservoirs. 
At low energies each reservoir constitutes one screening 
channel for the local spin thus leading to an over-screened 
situation for K > 1. We will thus use the terms reservoir 
and channel interchangeably. Furthermore, the spin-1/2 
is subject to an external magnetic field /ig. Each reser- 
voir consists of a left (L) and right {R) lead which are 
held at chemical potentials /xj^,, a — L,R, i — 1,. . . ,K 
(see Fig. [T]). Specifically we consider the system 

ii = ^ efe 4^,,^c,aka + hQ 

ioika , ^ 

Here c]^^^^ and Ciaka create and annihilate electrons with 
momentum k and spin a =t, i in lead a of the reservoir 
«, a denotes the Pauli matrices, and S is the spin-1/2 
operator on the dot. The antiferromagnetic exchange 
coupling Jo > is dimensionless in our convention. We 
stress that the exchange term does not couple differ- 




FIG. 1. (Color online) Sketch of the if -channel Kondo model 
([T]). The channels are labeled by i — I, . . . ,K, each channel 
is decomposed into two leads a = L,R held at chemical po- 
tentials l-ii^/n = ±Vi/2. The spin-1/2 on the dot is subject 
to a magnetic field ho. The exchange interaction J between 
dot and leads present in all channels is depicted for the Kth 
channel explicitly. 



ent reservoirs. The chemical potentials in the leads are 
parametrized by M^/^ = thus applying a bias volt- 

age Vi to each reservoir. Furthermore we introduce an 
ultra-violet cutoff D in each reservoir via the density of 
states N{uj) = vqD^ /{D"^ +uP). In the absence of a mag- 
netic field the transport properties of the model ([T]) have 
been previously studied in Ref . IHl 

RG analysis. — Following Ref. HU we have performed a 
two-loop RTRG analysis including a consistent deriva- 
tion of the relevant relaxation rates. This method has 
been successfully applied to study transport prope rties 
of other Kondo-type quantum dots in the past .1^211] xhe 
starting point is an RG equation for the renormalized 
exchange coupling J(A) obtained by integrating out the 
high-energy degrees of freedom in the reservoirs. To ac- 
complish this the one introduces a cutoff A into the Fermi 
function and integrates out the Matsubara poles on the 
imaginary axis by decreasing the cutoff from its initial 
value Ko ^ D down to some physical energy scale. The 
RG equation for the _ftr-channel model ([l]) reads up to 
two loop 



dA 



I3{J) = -2J2(1 - if J), 



(2) 



which defines the reference solution J(A) for our analysis. 
The RG flow has the well-knowrP non-trivial fixed point 
J* — 1/K. The scaling dimension of the leading irrele- 
vant operator is A = P'{J*) = 2/K, which is valid for 
K ^ 1 while the exact result is given 
The RG equation possesses the invariant 
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FIG. 2. (Color online) Differential conductance of an individ- 
ual channel for identical bias voltages Vi = V and K — 10 
(dashed lines) or if = 20 (solid lines) channels. We use the 
normalization Ga ~ (^)^^A^. For zero field one find^ 
G/Ga = 1 - 2{V/Tk)'^, I.e. G/Ga ^ 1 for V 0. Inset: 
Comparison of RTRG (solid line) and renormalized perturba- 
tion theory (PT) up to 0{,Jl) (dashed line) for = Tk and 
K — 20. We note that the broadening of the feature a.tV = h 
is not captured in PT. 



renormalized exchange coupling Jc = J(Ac) at the phys- 
ical energy scale^ (we consider T — 0) 



A, = : 



c{V,ho}. 



(6) 



Tk = An 



eJn 



J*- Jo 



K/2 



,-l/2Jo 



(3) 



which defines the Kondo temperature. With the initial 
condition Jq = J{Aq) the solution of the RG equation 
can be explicitly written as 



J(A) 



J* 



1 + Wiz)' 



-V 



(4) 



Here W{z) denotes the Lambert W functiod^ de- 
fined by z = W{z)e^''^\ which satisfies W'{z) = 
W{z)/z/{l + W{z)) for z ^ 0. The fixed point J* is 
reached for A —> as W{0) = 0. We note that the solu- 
tion ([4| is valid for all Jq but, of course, the derivation 
of ([2]) requires Jo <C 1. If Jq < J* we can perform the 
scaling limit Ag — > 00 and Jq — > while keeping the 
Kondo temperature constant. In this limit the solution 
for A <C Tk simplifies to 



J(A) = J* 



1 - 



A 



(5) 



i.e. there is a characteristic power-law behavior. Observ- 
ables are calculated in a systematic expansion around the 
reference solution and can be expressed in terms of the 



In contrast to the one-channel Kondo model the exis- 
tence of the attractive fixed point J* = 1/K implies that 
this expansion is well-defined for all Ac provided K ^ 1. 
We have calculated the effective dot Liouvillian and the 
current kernel yielding the renormalized magnetic field, 
the spin relaxation rates, the dot magnetization, and the 
current including the leading logarithmic corrections. All 
calculations follow Ref. llj we present here only the re- 
sults and discuss their properties. 

Identical bias voltages. — First let us assume that all 
bias voltages are identical, i.e. Vi = V. Straightforward 
calculation yields the renormalized magnetic field 



h=[l-K(J,-Jo)]h„, 



(7) 



which gives for the g factor in leading order (we consider 
the scaling limit Jq ^ from now) 



A 



(8) 



The longitudinal and transverse spin relaxation rates 



Fi =^ 



F2 = 



h+^{\V -hl^ + V + h) 



V + h+-{\V-h\^ + V + h) 



KJl (9) 
KJl (10) 



2 



with 



1,2 X 
\x\, — — xarctan — . 
TT Ti ' 



1 = 1,2. 



(11) 



Wc note that for Ac ^ Tk the renormalization of the 
magnetic field Q is much stronger than in the one- 
channel Kondo model, as the spin on the dot is coupled 
to more screening channels. 

Explicit formulas for the dot magnetization and the 



current are given in Eq. ( 14 ) below for the case of dif- 
ferent bias voltages. Simplifying io Vi = V we obtain 
the differential conductance Gi = dli/dVi per channel, 
which is plotted in Fig. [2] The conductance sharply in- 
creases around V = h where inelastic cotunneling pro- 
cesses start to contribute. Due to the strong renormaliza- 
tion of the magnetic field the increase occurs at voltages 
much smaller than the applied magnetic field Hq. Close 
to the resonance we find 



'lJ'^[l + 2J^C2{V -h)] iorV<h, 
^t:JI[1 + 2ttJcC2{V -h)] forV^>/i, 



where we assumed \V — h\ <^ h and defined 

A, 



CAx) = In 



1 = 1,2. 



(12) 



(13) 



We note that (12 1 is formally identical to the differen- 
tial conductance in the one-channel Kondo modeP^^ ex- 
cept for the functional form of the renormalized coupling 
([5]). In particular, there is no explicit dependence on the 
number of channels. We further note that close to the 



fixed point the conductance (12 1 is a quantity of order 
The logarithmic divergencies aX V = h are cut 
off by the transverse spin relaxation rate Thus, as 
r2 ~ ^/K, this feature becomes sharper with increasing 
K. At small voltages V < h the conductance is inde- 
pendently of the bias voltage given by Gi = nJ^/A. In 
particular, a power-law behavior in the voltage is only 
found for vanishing magnetic field^ while ([s]) yields a 
power-law dependence on the magnetic field. In the 
limit of a large field /i ^ Tk the linear conductance 
can be derived usin^i^ W{z) In z for z — > 00 to be 
Gi{V = 0,/i) = 7r/16/ln^(/i/rK), which is identical to 
the result in the one-channel model. 

Different bias voltages. — In the following we relax the 
condition Vi = V and consider channel-dependent bias 
voltages Vi. This introduces further energy scales which 
will give rise to additional features in the dot magnetiza- 
tion and thus the differential conductance. In this gen- 
eral set-up the magnetization and current through the 
ith channel are given by 



M 



A 
2/2' 



= + 2Af /, 



M 



(14) 



where /i_2 denote rates appearing in the Liouvillian and 
/I are similar terms in the current kernel. Explicitly 
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FIG. 3. (Color online) Differential conductance d through 
channel i for ho = 100 Tx and K = 20. The bias voltages 
are applied such that series of five channels have identical 
voltages, i.e. Vi = aiV with ai,...,5 = 1, as,.. .,10 = 1/2, 
an,. ...15 = 1/3, and ai6,...,20 = 1/4. For ho = the con- 
ductance is similar to the black curve in Fig. [2] i.e. there are 
no features. Inset: Same situation as in the main panel, but 
plotted with respect to the respective bias voltage Vi. 



these rates read 
/i = 27rKJ^h + AT:KJlhC2{h) 

-27rJf E,.i...k(^^ - h)C2{V, - h), (15) 
/2 = nKJ^h + 2Tr KJ^hC2{h) 

\jcT.,^,,„i,i\^^-''\2 + V^+h) (16) 

E,.i...K ( 1^^- - ^^12 - + - h), 

f} = \'^JlV,+TTJlV,Ci{V,) 

+TTjl[V,-h)C2{V,-h), (17) 

/m - -\Jl{ \v. - h\^ -V,-h)+ T^J^V^C, {v^ 

+^JlhC2{h) - \V, - h\^C2{V, - h). (18) 

The relaxation rates for the case of different bias volt- 
ages are obtained by straightforward generalization of 
Eqs. (|9| and ([lO|. In the derivation of Eqs. ([T5|-([l8]) 
we have neglected all terms in order that do not con- 
tain logarithms at either h = V = ox V = h. Thus 
when calculating the magnetization one has to expand 
consistently up to this order. 

We stress that although the different channels are 
not directly coupled via exchange interactions, the fi- 
nite magnetization M mediates a feedback between them. 
Consider for example the differential conductance of 
channel 1, i.e. the black curve in Fig. [U The sharp in- 
crease dXV = Vi = his again due to the onset of inelastic 
cotunneling processes. However, there is a second feature 
aXV~2h which is caused by the non-trivial voltage de- 
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FIG. 4. Derivative of the magnetization with respect to the 
voltage, i.e. dM/dV. All parameters are as in Fig. |3] We 
observe pronounced features a.t V = h and V — 2h. Inset: 
Magnetization as a function of the bias voltage. 



pendence of the dot magnetization. This can be seen in 
the derivative dM/dV shown in Fig. |4] which directly 
enters the differential conductance [see ([m])]. Similarly, 
the conductance in channel 6 (green curve in Fig. [s]) pos- 
sesses a feature dX V — 2Vq = h which is caused by the 
effect of the applied voltage in channels 1 to 5 onto M, 
while the increase a,t Vq = h {i.e. V — 2h) is due to the 
onset of inelastic cotunneling in channel 6. In this way 
the nonequilibrium magnetization introduces additional 
features into the individual conductances. We stress that 
such coupling effects between the channels are absent for 
vanishing magnetic fielcP or if all applied bias voltages 
are identical (see Fig. [2]). 

As a special case of the general set-up with channel- 
dependent bias voltages we can recover the experimen- 
tal situation realized- by Potok et al. in a semicon- 



ductor quantum dot. This is achieved by setting the 
chemical potentials to = ±F/2 and /i^/^ = for 

i — 2, . . . , K . In each of the channels 2, . . . , K we intro- 
duce even and odd combinations of the electron opera- 
tors, {ciLkcr^Cifika) / V2, such that the even combinations 
couple to the spin on the dot with exchange interaction 
2J while the odd ones decouple completely. The resulting 
RG equation for J is given by ^ and thus possesses the 
fixed point J* . The experimental set-up^ is now obtained 
by specializing to = 2. In the presence of a magnetic 
field the resulting differential conductance is very similar 
to the case of identical bias voltages shown in Fig. [2j In 
particular, since the second channel does not provide an 
additional energy scale there appear no features in the 
differential conductance beside the cotunneling peak at 
V = h. A possible experimental set-up to observe the 
additional features shown in Fig. |3] thus requires at least 
two channels with non-zero and different bias voltages. 

Conclusions. — To sum up, we have studied the 
nonequilibrium transport properties of a multichannel 
Kondo quantum dot in the presence of a magnetic field. 
We used the solution of the two-loop RG equation to 
derive analytical results for the g factor, the spin relax- 
ation rates, the dot magnetization, and the differential 
conductance. The latter shows typical features of inelas- 
tic cotunneling. We showed that the main difference to 
the previously studied^ situation without magnetic field 
is the appearance of additional features in the differential 
conductance, which originate in the feedback between the 
channels mediated by the finite dot magnetization. 
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